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ABSTRACT 

Previous clinical studies in pulmonary arterial hypertension (PAH) have concentrated predominantly on distal pulmonary vascular 
resistance, its contribution to the disease process, and response to therapy. However, it is well known that biomechanical factors 
such as shear stress have an impact on endothelial health and dysfunction in other parts of the vasculature. This study tested 
the hypothesis that wall shear stress is reduced in the proximal pulmonary arteries of PAH patients with the belief that reduced 
shear stress may contribute to pulmonary endothelial cell dysfunction and as a result, PAH progression. A combined MRI and 
computational fluid dynamics (CFD) approach was used to construct subject-specific pulmonary artery models and quantify 
flow features and wall shear stress (WSS) in five PAH patients with moderate-to-severe disease and five age- and sex-matched 
controls. Three-dimensional model reconstruction showed PAH patients have significantly larger main, right, and left pulmonary 
artery diameters (3.5 ± 0.4 vs. 2.7 ± 0. 1 cm, P = 0.01 ; 2.5 ± 0.4 vs. 1 .9 ± 0.2 cm, P = 0.04; and 2.6 ± 0.4 vs. 2.0 ± 0.2 cm, P = 0.01 , 
respectively), and lower cardiac output (3.7 ± 1 .2 vs. 5.8 ± 0.6 L/min, P = 0.02.). CFD showed significantly lower time-averaged 
central pulmonary artery WSS in PAH patients compared to controls (4.3 ± 2.8 vs. 20.5 ± 4.0 dynes/cm 2 , P = 0.0004). Distal 
WSS was not significantly different. A novel method of measuring WSS was utilized to demonstrate for the first time that WSS is 
altered in some patients with PAH. Using computational modeling in patient-specific models, WSS was found to be significantly 
lower in the proximal pulmonary arteries of PAH patients compared to controls. Reduced WSS in proximal pulmonary arteries 
may play a role in the pathogenesis and progression of PAH. This data may serve as a basis for future in vitro studies of, for 
example, effects of WSS on gene expression. 

Keywords: biomedical engineering, computer simulation, endothelium-derived factors 



Pulmonary arterial hypertension (PAH) is a chronic disease 
characterized by progressive elevation in pulmonary artery 
pressure due to increased resistance in the distal pulmonary 
arteries and decreased compliance in the proximal arteries. 
If left untreated, PAH patients have a median survival of 
less than three years from the time of diagnosis, 1121 often 
attributable to the resulting compromised right ventricular 
function. A known mechanism of endothelial dysfunction 
in the systemic circulation is reduction in wall shear 
stress [WSS]; however, little is known about pulmonary 
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endothelial cell-specific responses to these forces and their 
role in the pathogenesis and progression of pulmonary 
vascular diseases, such as PAH. [341 It has been proposed 
that pulmonary vascular endothelial dysfunction may play 
a large role in the pathogenesis of PAH, [5 91 including the 
medial thickening, intimal fibrosis, microthrombosis, and 
plexiform lesions that affect the pre-capillary pulmonary 
arteries. 
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Previous clinical studies have concentrated on the properties 
of the distal vasculature in PAH, its contribution to the 
disease process, and its response to therapy. Specifically, 
pulmonary vascular resistance is routinely measured 
and used as a guide for therapeutic efficacy and therapy 
escalation. 110111 Ventricular loading, however, is dependent 
on both resistance and total arterial compliance. [1214] 
Decreased compliance associated with PAH and changes in 
response to vasodilators have been demonstrated in patients 
with PAH. [15161 Despite this and other compelling evidence 
that pulmonary artery compliance is a significant factor in 
decoupling the right ventricle from its vascular load, [1718] few 
clinical investigations to date have looked atbiomechanical 
factors, such as wall shear stress, in the central pulmonary 
arteries in PAH, [19] and none have examined their potential 
to offer insight into disease severity and outcomes. 

The ability to accurately measure WSS through the 
combination of MRI and computational fluid dynamics 
(CFD] is a novel development. WSS, measured in dynes/cm 2 , 
is a frictional force along the inner wall of the artery caused 
by viscous drag. Unlike ultrasound, which only measures 
one component of velocity, CFD allows for the creation 
of a time-varying 3D velocity field within the lumen. CFD 
facilitates analysis of flow patterns, spatial distribution 
of velocity, turbulence, and WSS in a diseased state. To 
improve accuracy in the setting of complex geometry, 
"meshing" is used, in which the arterial geometry is broken 
down into smaller components, such as tetrahedrons. 1201 
High-resolution 3D imaging, such as gadolinium-enhanced 
magnetic resonance angiography, combined with CFD, 
enable accurate assessment of WSS in the pulmonary 
vasculature as previously demonstrated. [21] 

We hypothesize that the abnormal anatomic and 
hemodynamic properties of the central pulmonary arteries 
will lead to a decrease in WSS which, as demonstrated in 
the systemic circulation, may lead to pulmonary endothelial 
cell dysfunction and abnormal gene expression. As a first 
step toward elucidating the role of shear stress in the 
pathogenesis and progression of PAH, we combined MRI and 
computational fluid dynamics to construct subject-specific 
pulmonary artery models in five PAH patients and five 
age- and sex-matched controls in order to quantify WSS 
and provide a basic range of normal WSS for future studies. 

MATERIALS AND METHODS 

Magnetic resonance imaging 

Gadolinium-enhanced magnetic resonance angiography 
(MRA) scans of five subjects with PAH (16- and 41-year-old 
males, and 19-, 36- and 50-year-old females] and five 
age- and sex-matched healthy subjects in the supine 
position were obtained in a 1.5T GE Signa MR scanner 



(GE Medical Systems, Milwaukee, Wise.]. IRB-approved 
informed consent was obtained from all subjects. 
Gadopentetate dimeglumine contrast agent (Magnevist, 
Bayer HealthCare, Wayne, N.J.] 0.2 mmol/kg was injected 
at a rate of 2 mL/sec via antecubital intravenous access, 
and a three-dimensional fast gradient echo MRA sequence 
was used to obtain a volume of coronal slices in 20 seconds 
during one breath-hold. Slice thickness ranged from 3.0 
to 3.2 mm, interpolated to 1.5 to 1.6 mm, totaling 50 to 
60 slices per volume. A 512 x 192 acquisition matrix 
(reconstructed to 512 x 512]was used with an in-plane field 
of view of 3 5 x 3 5 cm 2 to provide an in-plane reconstructed 
spatial resolution of 0.68 mm. The true imaging resolution 
of the MRA scans was 0.7 x 1.8 x 3.0 mm while the 
interpolated image voxel size was 0.7 x 0.7 x 1.5 mm. Other 
scan parameters included repetition times (TR] of 2.7 to 
3.5 milliseconds, echo times (TE] of 0.704 to 0.844 ms, and 
flip angles of 25°. 

During the same imaging study, breath-held, segmented 
k-space phase contrast cine images with through-plane 
velocity encoding were obtained perpendicular to the 
main (MPA], left (LPA], and right (RPA] pulmonary arteries. 
This series of images provided subject-specific physiologic 
flow information at the three locations. The acquisitions 
were ECG-gated, with TR = 2.2-2.5 ms and eight views 
per segment, resulting in a temporal resolution of 40 
ms. Twenty time points evenly spaced within a cardiac 
cycle were retrospectively reconstructed. Through-plane 
velocities were measured through 10 -mm thick slices with 
velocity encoding gradients of 1 5 0 cm/ seconds to 2 5 0 cm/ 
seconds. Other imaging parameters included 28 or 30 
square cm field of views, a 256 x 160 acquisition matrix 
yielding an in-plane resolution of 1.1 x 1.8 mm, 4.4 to 4.7 
ms TRs, 2.2 to 2.5 ms TEs, and 15° or 20° flip angles. 

Image processing, model construction and 
boundary condition specification 

In each volume of MRA data acquired, geometric distortions 
caused by magnetic gradient nonlinearity were corrected 
to preserve the true dimensions of the pulmonary 
vasculature. 1221 Three-dimensional, subject-specific solid 
models of the pulmonary tree were created from the MRA 
images using custom software 1231 and discretized into a 
coarse, isotropic mesh using a commercially available, 
automatic mesh generation program (MeshSim, Simmetrix, 
Clifton Park, N.Y.). 

The cine phase contrast images (PC-MRI] taken 
perpendicularly to the main, left, and right pulmonary 
arteries were used to calculate in vivo, time-resolved 
volumetric flow. For each of the 20 time points, the 
lumen of the respective vessels was determined from the 
intensity magnitude images using a level set method, 1241 
and through-plane velocity values for each pixel were 



471 



I 



Tang et al.: Shear stress in PH 



integrated over the bounded area to calculate a total volume 
flow rate. 1221 The periodic flow waveform at the MPA was 
mapped onto the inlet face using a parabolic profile, and 
the average volumetric flow calculated from the PC-MRI for 
the left and right pulmonary arteries was used to determine 
the flow split to each side. Resistance values were assigned 
at the outlet of each branch vessel based on the area of 
the outlet, the LPA/RPA flow split, and either healthy 
pressure values (average of 12 mmHg] for the healthy 
subjects or subject-specific pressure values measured from 
catheterization in the PAH subjects. Furthermore, a rigid 
wall assumption and a no-slip boundary condition at the 
wall were also prescribed. 

Once all boundary conditions were applied, the flow 
solution for each subject-specific mesh was obtained 
using a stabilized finite element method to solve the 
incompressible Navier-Stokes equations. 1251 Each coarse, 
isotropic mesh was then adaptively refined into an 
anisotropic mesh based on the results of the solution (in 
order to reduce local error), and this process was repeated 
until mesh independence was reached. The resulting 
refined subject-specific meshes contained an average size 
of 1.2 million tetrahedral elements. 

Analysis of simulation results 

From the simulation results, time-averaged wall shear 
stress (x ] was computed for each subject, t is defined 

^ mean-' ~ ' mean 

as the magnitude of the time-averaged surface traction 
vector, t^ which is the tangential component of the traction 
vector, t. 




where: 

t. = t- (t« n]n 
and: 



t = an 

where s is the stress tensor, n is the surface normal vector, 
and T is the period of the cardiac cycle. 

In order to quantitatively analyze and compare the mean 
wall shear stress for all subjects, 10-mm strips were taken 
perpendicularly around the circumference of the LPA 
and RPA before the first branch, as well as at a random 
selection of distal branch vessels (three vessels on each 
side with cross-sectional area of 10-mm 2 ], and values were 
averaged over the surface area. Mean wall shear stress 
was statistically compared using a student's t-test with 
significance set at P < 0.05. 

RESULTS 

Anatomic and hemodynamic features of PAH 

Three-dimensional model construction of the large 
pulmonary arteries and major branches illustrated the 
anatomic differences between healthy, normal subjects, and 
PAH patients. As predicted, PAH patients were observed 
to have more tortuous branch vessels and larger proximal 
arteries, with an average main pulmonary artery diameter 
of 3.5 ± 0.4 cm (compared to 2.7 + 0.1 cm in healthy subjects, 
P = 0.01], right pulmonary artery diameter of 2.5 ± 0.4 cm 
(compared to 1.9 ± 0.2 cm in healthy subjects, P = 0.04], and 
left pulmonary artery diameter of 2.6 ± 0.4 cm (compared to 
2.0 ± 0.2 cm in healthy subjects, P= 0.01]. Additional subject 
data collected during the study is summarized in Table 1. 

Flow rate (i.e., cardiac output] measured at the MPA was 
significantly lower in the PAH patients (3.7 ± 1.2 L/min 
(2.0 L/m/m 2 ] compared to 5.8 ± 0.6 L/min (3.3 L/m/m 2 
in healthy subjects, P = 0.02]. The combination of larger 
proximal arteries and lower flow rates in the PAH patients 
led to lower velocity magnitude and more stagnant flow 
throughout the cardiac cycle in the proximal arteries 
of the PAH patients. A volume-rendered plot of velocity 



Table 1 : Subject-specific MRI and catheterization data 



Subject Age Sex BSA MRI data Cardiac catheterization data 

( m2 ) MPA diameter MPA flow RA p (mean, PA P (S/D, M; PVRI 

(cm) (L/min) mm Hg) mm Hg) (U/m 2 ) 

PAH 1 16 M 1.55 3.2 2.6 4 70/37, 50 10.7 

patients 2 19 F 1.58 3.1 2.4 7 95/45,67 21.5 

3 36 F 2.16 3.2 5.2 4 87/36, 50 22.9 

4 41 M 1.87 3.7 4.4 4 77/30, 49 18.1 

5 50 F 2.00 4.1 3.8 7 110/42,68 37.2 
Controls 21 M 1.76 2.7 6.0 

22 F 1.66 2.7 5.1 

31 F 1.60 2.6 5.4 

35 M 2.00 2.6 6.6 

51 F 1.75 2.8 5.9 



MPA: main pulmonary artery; RA p : right atrial pressure; PA p : Pulmonary artery pressure (systolic/diastolic, mean); PVR: pulmonary vascular resistance 
(Woods units) 
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magnitude in the pulmonary arteries of one representative 
PAH patient and an age- and sex-matched normal subject 
is shown in Figure 1. Note that, even at peak systole, 
the velocity magnitude in the proximal pulmonary 
arteries of the PAH patient was very low, and almost no 
flow movement was observed throughout diastole. In 
comparison, the velocity magnitude was consistently 
higher in the pulmonary arteries of the normal control 
subject at similar time points during the cardiac cycle. 
A more detailed look at the cross-sectional flow features in 
the LPA and RPA of these two subjects is shown in Figure 2. 
A low-intensity, more uniform flow field was observed in 
the PAH patient whereas high-intensity, swirling flow along 
the walls of the proximal arteries was found in the normal 
control subject. 

Mean wall shear stress (WSS] results are summarized for all 
subjects in Figure 3. Plots of mean WSS illustrate that PAH 
subjects have lower mean WSS in their proximal arteries 
compared to healthy subjects. In many subjects, areas of low 
mean WSS were found at sites opposite to branch points, 
indicating possible recirculation zones in these regions. 
However, it should be noted that the levels of low mean 
WSS in the proximal pulmonary arteries of the control 
subjects is still higher than the overall levels found in the 
PAH patients. When quantified using 10-mm strips taken 
around the circumference, the mean WSS of the proximal 
pulmonary arteries of healthy subjects was nearly five times 
higher than that of the PAH patients (20.5 ± 4.0 dynes/cm 2 
in the control subjects compared to 4.3 ± 2.8 dynes/cm 2 



in the PAH patients, P = 0.0004]. Mean WSS values in the 
distal arteries was 1.4 times higher in the normal subjects 
than the PAH patients (14.1 ± 0.7 dynes/cm 2 in the control 
subjects compared to 10.1 ± 0.9 dynes/cm 2 in the PAH 
patients, P= 0.004]. 

DISCUSSION 

Altered systemic shear stress has been shown to have 
a detrimental effect on cardiovascular health. As an 
initial step in elucidating whether similar biomechanical 
forces adversely affect pulmonary vascular health, 
we quantified the mean WSS in patients with PAH as 
compared to age- and sex-matched controls. We found 
significantly lower wall shear stress in the pulmonary 
arteries of five subjects with moderate or severe PAH 
compared to control subjects using combined magnetic 
resonance imaging and computational fluid dynamics to 
visualize and quantify subject-specific, three-dimensional 
hemodynamic conditions. This study also provides novel 
insight into the time-varying blood flow features that exist 
in the pulmonary vasculature in health and disease and 
suggests a mechanism by which conditions of low shear 
stress may result in disease progression. 

Numerous secondary observations can be made about 
alterations in velocity and flow in the pulmonary 
vasculature of PAH. The swirling, high-wall-velocity flow 
noted in normal subjects at rest did not occur in PAH 



Pulmonary Arterial Hypertension Patient 




0 cm/sec Velocity Magnitude 30 cm/sec 



Figure 1 : Volume-rendered velocity magnitudes plotted in the pulmonary arteries of a 50-year-old female with Pulmonary Arterial Hypertension (top) and a 
5 1 -year-old healthy female volunteer (bottom). The left-most panel shows peak-systole for both subjects, and equally-spaced time points throughout the cardiac 
cycle are shown in the following three panels. 
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Figure 2: Cross-sectional view of flow features in the LPA and RPA of the 
PAH patient (right) and the age- and sex-matched normal control subject 
(left) featured in Figure 1. Arrowheads (shown enlarged in inset of panel b) 
are used to represent velocity vectors and are colored and scaled according 
to velocity magnitude. In addition, velocity magnitude is plotted on a cross 
sectional slice through the vessel. 

patients, thus leading to significantly lower wall shear 
stresses. The most notable difference in flow features 
between a PAH patient and a normal subject, as shown in 
Figures 2 and 3, are the higher, swirling velocities along the 
wall in the normal subject. In comparison, flow is basically- 
stagnant for more than half of the cardiac cycle in the PAH 
patient, with higher flow velocities along the wall only 
occurring at peak systole. These higher flow velocities along 
the walls of the proximal vessels in normal subjects led to 
an overall six times higher wall shear stress than in PAH 
patients when compared across all subjects. Interestingly, 
this difference in wall shear stress was diminished in the 
distal vasculature, likely due to pruning of the distal arteries, 
thus leading to similar distributions of flow. 



It has been extensively demonstrated in the systemic 
circulation that wall shear stress modulates endothelial 
function. Specifically, systemic increases in WSS result 
in elevated expression of eNOS. [3;26_28] Nitric oxide is also 
a primary modulator of pulmonary vasodilation and 
target of PAH drug therapies, such as sildenafil. [2932] In 
addition, WSS alters pulmonary expression of Hsp90, an 
activator of eNOS, NAD(P)H dehydrogenase quinone 1, an 
antioxidant, monocyte-chemotactic protein 1, and platelet 
endothelial cell adhesion molecule 1.P 8 ' 33 - 36 ] Reduction 
in angiotensin-converting enzyme has also been shown 
with increases in pulmonary vasculature shear stress. [37] 
Therefore, one could conclude that dramatic changes in 
the WSS of the pulmonary vasculature, as we demonstrated 
in PAH, would have significant effects on pulmonary 
endothelial health and remodeling. 

Our findings may elucidate one of the mechanisms by 
which epoprostenol improves patient outcomes, results 
in long-term reductions of pulmonary vascular resistance, 
and partially reverses disease progression. 1381 Epoprostenol 
has been shown to significantly affect cardiac output and 
pulmonary vascular resistance and to cause long-term 
vasodilation despite lack of response to acute vasodilator 
testing. This may suggest a unique mechanism of action 
when given chronically. 1381 We propose epoprostenol 
may alter wall shear stress in the pulmonary arteries by 
increasing cardiac output resulting in a positive remodeling 
response and long-term vasodilation. However, our current 
therapies at best stabilize but do not reverse the process. 
Perhaps this frustration is because we are combating two 
entities: the disease itself and the resulting abnormal 
biomechanical effects on the endothelium and gene 
expression. If shear stress is found to be a major determinant 
of and contributor to PAH, a target of future therapies may 
be alteration of the biomechanical environment as a means 
of reversing disease progression. Great interest now exists 
in using more aggressive therapy, e.g., epoprostenol, earlier 
in the disease treatment algorithm to get the most effect 
possible early on. This study may lend additional support 
to this proposed approach. 

Clinically, assessment of wall shear stress by magnetic 
resonance, an imaging modality with increased utilization 
in clinical practice for the evaluation of right ventricular 
function, may serve to stratify disease severity and monitor 
the efficacy of existing and emerging therapies. Currently, 
symptomatology and PVR are the primary methods 
of assessing disease severity and response to therapy. 
However, measurement of WSS may prove to be a more 
comprehensive, non-invasive metric for tailoring therapy 
and determining prognosis. 

The primary limitations of this study are the small number 
of patients and the severity of their disease. While there was 
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Figure 3: Summary of time-averaged wall shear stress (WSS) results for all subjects. Mean WSS color maps (left panel) for two representative age- and 
sex-matched (top) normal subjects and (bottom) PAH patients. Mean WSS was lower in the proximal arteries of the PAH patients than in those of the normal 
subjects. Mean WSS averaged over the area of 10-mm circumferential strips taken at the LPA and RPA (proximal) and distal locations was significantly different 
between the two populations (right panel). 



an order of magnitude difference in WSS for each patient 
compared to control, this represents a small study of only 
five severe PAH patients. Of note, the broad age range of 
patients (16 to 50) was purposely chosen to increase the 
heterogeneity of the patient population and, therefore, 
improve possible applicability of findings. While this pilot 
study found a significant difference in WSS in PAH patients, 
and may lead to further advances in the basic science of 
this disease, future work must investigate the use of these 
non-invasive imaging and CFD techniques in a greater 
number of patients including a "healthier" PAH population 
to determine the sensitivity and specificity of these methods 
and clinical applicability. 

This study represents the first time that three-dimensional 
hemodynamic conditions have been quantified in 
patient-specific models of pulmonary arterial hypertension 
and is an important step toward elucidating the role of 
shear stress in regulating pulmonary artery endothelial cell 
dysfunction. Our results point to a biomechanical mechanism 
as a component of pulmonary endothelial cell dysfunction 
and further progression of the disease. As distal resistance 
increases, flow is decreased in the enlarged proximal vessels, 
leading to decreased wall shear stress, which is known to 
negatively impact endothelial cell health. This shear-mediated 
endothelial dysfunction can lead to further remodeling of 
the proximal vessels, as well as downstream remodeling via 
endocrine signaling, resulting in worsening of the disease and 



further decreases in wall shear stress. All of these changes 
suggest a contribution of the unhealthy disease state to 
disease progression, that is, PH begets PH. 
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